Radio frequency spectroscopy of a strongly imbalanced Feshbach-resonant Fermi gas 
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A sufficiently large species imbalance (polarization) in a two-component Feshbach resonant Fermi 
gas is known to drive the system into its normal state. We show that the resulting strongly- 
interacting state is a conventional Fermi liquid, that is, however, strongly renormalized by pairing 
fluctuations. Using a controlled 1/N expansion, we calculate the properties of this state with a 
particular emphasis on the atomic spectral function, the momentum distribution functions displaying 
the Migdal discontinuity, and the radio frequency (RF) spectrum. We discuss the latter in the light 
of the recent experiments of Schunck et al. (Science 316, 867 (2007)) on such a resonant Fermi 
gas, and show that the observations are consistent with a conventional, but strongly renormalized 
Fermi-liquid picture. 

PACS numbers: 67.85.De, 03.75.Kk, 03.75.Ss 



I. INTRODUCTION 

One of the key recent developments in studies of degen- 
erate atomic gases is the tunability of atomic interactions 
via a Feshbach resonance (FR). This has led to the real- 
ization of a resonantly-paired s-wave superfluid that can 
be tuned between the two well-studied limits of a weakly- 
paired Bardeen-Cooper-Schrieffer (BCS) superfluid and 
a strongly-paired diatomic molecular Bose-Einstein con- 
densate (BEG) superfluid 1, 2, 3, 4,i,[|. 

The two asymptotic superfluid regimes at large posi- 
tive and negative FR detunings allow a detailed quan- 
titative description, made possible by the existence of a 
small gas parameter, na^, corresponding to the ratio of 
a short scattering length a to a large average particle 
spacing i = n^^/'^. In contrast, although the intermedi- 
ate low temperature crossover regime is a conventional 
superfluid, that smoothly interpolates between the BGS 
and BEG limits, its quantitative description (in a broad 
resonance case) is hindered by strong interactions, char- 
acterized by a diverging scattering length and absence of 
a natural small parameter Q. The flip side, of course, is 
that a diverging scattering length leaves particle spacing 
as the only relevant length scale, leading to a universal 
phenomenology of a resonant Fermi gas near a unitary 
point. 

While earlier studies focused on the case where the 
populations of the two atomic species involved in pair- 
ing are equal (vanishing polarization), and thus on the 
nature of the superfluid phase, recent experimental and 
theoretical investigations have extensively explored the 
imbalanced resonant Fermi gas, extending its FR detun- 
ing phase diagram to a finite polarization [HI [13, [a, 
Il4l . [isl . [igI . [rn | . These studies have consistently found 
that at high polarization the state is non-superfluid, and 
have treated it as a simple Fermi gas. However, at low 
detuning the system is strongly resonantly interacting. 



and a detailed description of the highly polarized normal 
state near the resonance remains a challenge [l8l |. 

Some light on the complex nature of the strongly- 
interacting normal state has been shed by a recent radio 
frequency (RF) spectroscopy experiment [13] ■ At high 
polarization an absence of a BEG peak on one hand and 
presence of a temperature-dependent spectral shift rela- 
tive to the atomic line, on the other, was observed. In- 
terpreting the latter as a pairing gap i , taken together, 
these observations have been interpreted as evidence for 
a paired non-superfluid state. While unsurprising at fi- 
nite temperature, and certainly present on the BEG side 
of the resonance, where the gap is set by molecular bind- 
ing energy, an existence of such a state at zero tempera- 
ture would constitute a dramatic departure from a mod- 
ern understanding of possible condensed matter ground 
states. In particular, the only established route to a sup- 
pression of boson (Gooper pairs or diatomic molecules 
in the present context) superfluidity at zero temperature 
is through a localization of bosons by quenched disorder 
into a Bose glass [l^l, by a commensurate (e.g., imposed 
optical) lattice into a Mott insulator [23, [2l[ , or by crys- 
tallization as in a case of a solid ''He. Since none of these 
mechanisms appears to be at play in the trapped dilute 
atomic gas studied in experiments by Schunck et al. [l9j , 
their observations and conclusions remain puzzling. 

The interpretation of RF spectra in the strongly- 
interacting regime has been the subject of numerous the- 
oretical investigations [13, [13, |m, 25, 2(1] , but remains 
incomplete. The recent RF spectroscopy experiments on 
a strongly polarized Fermi gas [13, [13], have rekindled 
theoretical studies of such a system with a focus on the 
unitarity regime 




29 . 13C 1 a nd its strong interactions 
"^'3|ll[33,[3i,[3i. 



in the normal state 



The purpose of the present paper is to make a case 
for a more conservative interpretation of the RF spectra 
of Refs. [13, [13, namely that a non-superfluid state of a 
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FIG. 1: Momentum distribution ncr{k) of the majority (|1)) 
and minority (|2)) atoms for polarization P = 0.9 at zero 
temperature and at resonance. The residues for the majority 
and minority atoms are, respectively, Zi = 0.56 and Z2 ~ 
0.29. 



highly polarized, resonantly interacting Fermi gas is in 
fact a Fermi liquid, albeit a strongly renormalized one. 
Indeed, the Luttinger relations discussed in Refs. [stI [ssI 
require any non-superfluid state to have Fermi surfaces 
for both the majority and minority species, enclosing the 
same volumes as for non-interacting fermions. We will 
show that such a Fermi liquid is perfectly consistent with 
the experiment of Ref. il9|, and that the observed shift 
in the RF spectrum of the minority atoms is the result 
of large self-energy effects due to the Feshbach resonant 
scattering, as also emphasized in parallel recent stud- 
ies [32,[3l|3l. Based on this we suggest that a far better 
test of the nature of the ground state is a measurement of 
the momentum distribution 712 (/c) of the minority atoms; 
in the following the label a = 1,2 denotes the majority 
and minority species, respectively. Such a measurement 
would be a direct test for the existence of a Fermi surface, 
marked by a Migdal discontinuity and characterized by 
the quasiparticle residues < < 1 



(1) 



illustrated in Fig. [TJ This hallmark Fermi liquid feature 
would be absent in the case of a paired ground state, 
allowing for a sharp qualitative distinction between two 
possibilities. An interesting feature of Fig. [1] is that the 
functions n„(k) are both universal as functions of k/kp 
at resonance i.e. for the Fermi gas at unitarity, the func- 
tional forms (including the values of the discontinuities 
Z„) depend only upon the polarization P. Our compu- 
tations of these universal functions is however not exact, 
and next we describe our computational method. 

We build on our earlier studies of the superfluid 
statejs^ H^: our theoretical approach to the treatment 
of this strongly interacting (small parameter-free) normal 
state is based on the introduction of an artificial small pa- 
rameter with N the number of distinct "spin"-l/2 
fermion flavors in the generalized model. The advantage 
of such a generalization is that, for N 00, the problem 
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FIG. 2: RF spectrum for polarization P = 0.97: the intensity 
13(1') (arbitrary units) vs. the detuning from the resonance 
frequency u measured in units of the Fermi energy ep . In free 
space, the resonance would be at = 0. The shift in the res- 
onance frequency above is due to strong interactions between 
fermions in the non-superfluid ground state of a polarized 
Fermi gas. Our primary claim is that such a shift is present 
even while the ground state remains a Fermi liquid, with the 
discontinuities in the momentum distribution function shown 
in Fig. m 



is exactly solvable, with finite N corrections computable 
via a systematic expansion in 1/N about this solvable 
limit. In particular, we calculate the RF spectrum to 
sublcading order in 1 /N using the full interaction matrix, 
including the appropriate vertex correction. We then 
use these controlled 1/N results to extrapolate to the 
experimentally-relevant case of a single flavor {N — 1) of 
two opposite- "spin" (hyperfine levels) fermionic atoms. 
This extrapolation is a subtle issue, and we cannot rule 
out the possibility of a nonanalyticity in the N —^ 1 limit. 
However, such nonanalyticities are rare and there is no 
reason to expect them here. 

We regard the large- TV expansion as providing a frame- 
work for understanding qualitative aspects of the renor- 
malized Fermi liquid properties of strongly-interacting 
imbalanced Fermi gases. A sample result for the RF spec- 
trum obtained via this approach is shown in Fig. [21 Note 
the asymmetric lineshape: this arises from the imagi- 
nary part of the pairing fluctuations propagator which 
contains the phase space for the decay of a Cooper pair 
into a two fermion final state. 

The universality considerations noted above for n„{k) 
also apply to the the RF absorption spectrum. However, 
here there are also additional complications associated 
with the scattering lengths of excited states, and these 
will be discussed below. 

The outline of the paper is as follows. In Sec. |TT] we 
introduce the resonant single-channel model as well as 
its A^-fiavor generalization, and carry out its systematic 
expansion in 1/A^. Then, in Sec. IIIIl we illustrate the 
Fermi liquid properties of the system at high polarization. 
In Sec. llVl the radio- frequency probe is described and we 
discuss our results in terms of the recent experiments of 
Schunck et ai, and conclude with a brief summary in 
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II. FORMALISM 



Model 



To capture the physics involved in the RF experiment, 
we model the system in terms of the three lowest hyper- 
fine states, |cr) = |f), |2} and |3). The two lowest states, 
|1) and |2), are loaded with atoms and responsible for 
the strong superfluid correlation, while the higher hyper- 
fine state, |3), is initially empty. Experimentally, the RF 
field at frequency uj is used to induce atomic transitions 
from the state |2) to the state |3), and the induced tran- 
sition rate is measured as a function of w. As we will 
show below, this allows one to experimentally probe a 
two-particle correlation function. 

Although our consideration for RF spectroscopy on 
fermionic atoms will be quite general, we will put par- 
ticular emphasis on recent experiments on ^Li. In this 
system, the three lowest lying states, |2) and |3) can 
be identified with \F = 1/2, = 1/2), \F = 1/2, mp = 
— 1/2) and \F = 3/2, mp = —3/2), respectively. 

By integrating out the higher order hyperfine states 
we obtain an effective Hamiltonian in terms of the three 
lowest states, given by 

3 

+ ^ E /rf'rA..,Vi(r)V'I.(r)V'..(r)V'.(r), (2) 

(T,cr' — 1 

where X^^i are couplings (interaction strengths set by the 
corresponding scattering lengths aaa-') between states \a) 
and \a'). ■!/'J(r), V'cr('') ^^^j respectively, the fermion cre- 
ation and annihilation operators at position r and hy- 
perfine state cr, which obey the usual anticommutation 

relation |^g.(r),'0^, (r')| = S{r ~ r')6cr.a' ■ The detuning 

of the level a is controlled by a Zeeman field encoded by 
the detuning parameter u)^. The chemical potential /in- 
fixes the average atom density in the hyperfine state 
(spin) a. 

While in above model the particle number in each hy- 
perfine state is a good quantum number, in principle 
there are additional interaction channels present in the 
physical system, that break this symmetry. For instance 
^ in 40js:, (where |1), |2), and |3) are \F = 9/2, mp = 
-9/2), |F = 9/2, mp = -7/2) and \F = 9/2, mp = 
—5/2)) a collision term such as '(/'i'03V'2V'2 conserves the 
total quantum number mp is therefore allowed by sym- 
metry. However, effects of such interactions, that do not 
conserve the atom number in each hyperfine state are en- 
ergetically suppressed by virtue of 2u)2 7^ -l-wa, due to 
the large quadratic Zeeman splitting between the three 



levels. Experimentally, the Feshbach resonances are cho- 
sen such that these kind of non-conserving processes are 
minimized, and therefore we will not consider them any 
further here. 

For a short-range s-wave interaction, the Pauli princi- 
ple enforces X^a = 0, which together with the exchange 
symmetry A^ct' — Xa-'a, reduces the nine coupling con- 
stants Acrcr' down to three. The corresponding three two- 
particle scattering lengths Uaa' for atoms in states \a) 
and \a') (for a ^ a') are related to the strengths of the 
couplings Acrcr' via the relation 



m 



1 



X„ 



d^k 1 



(3) 



where ek = fc^/(2m) is the the free fermion dispersion. 



B. Large TV expansion 

As discussed in the Introduction, because we are inter- 
ested in the system in the vicinity of the unitary point, 
where the atom interaction is strong, an analysis based on 
a straightforward perturbation theory in the gas param- 
eter na^ (that diverges at the unitary point) clearly fails 
and an approach nonperturbative in na^ is required. To 
this end, we employ the procedure, introduced and suc- 
cessfully utilized for the superfluid phase in Refs. [s^ HO, 
of generalizing the physical model to that of N flavors 
of each of the three hyperfine states species, with the 
Hamiltonian given by 

3 N 



3 N 

+ ^ E E y ^'rA..-^,L(r)vl,(r)V,.,,(r)V,.,(r). 

(7,a' — l — l 

(4) 

It is then straightforward to develop an expansion in 
powers of 1/A^, taking the physically relevant limit = 1 
at the end of the calculation. Since we are interested in 
the normal properties only, the full effective action for- 
malism of Refs. [3^ HO will not be required. Instead we 
may calculate the 1/N contributions diagrammatically, 
observing that vertices bring factors of 1 /N, and particle- 
particle loops a factor of N, as we sum over all compo- 
nents. Having thereby established the relevant class of 
leading 1/A'^ diagrams, for ease of notation, we then drop 
the i flavor index. 

In this way, we arrive at the thermal Green's function 
of species a defined by 



Ga{k, ILUn) 



dTe-"-(T,{V^,(k,r)V^t(k^O)}), 



(5) 

where w„ = 27r (n -|- i) //? is the fermionic Matsubara 
frequency and P — 1/T is the inverse temperature T and 
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FIG. 3: Self-energy diagram to leading order in 1/A''. 

Tt is the time ordering operator with the imaginary time 
T. As in the standard perturbation theory, we express 
the Green's function in term of its non-interacting form 
and the self-energy 

G;i(k,ic^„) = G^°)-i(k,ic^„) -E,(k,zc^„), (6) 

where So-lk, ?a;„) is the (Matsubara) self-energy, 

Gct"'' (k, ju;„) = (icj„ — ecr.k + /icr) ^ is the bare Green's 
function, and SaM = ek + tl^o- is the bare fermionic disper- 
sion. 



The self-energy is determined to lowest order in \/N . 
It can be expressed as 

3 

S^(k,iw„)= ^ (k, (7) 

a' = l(#<T) 

where is the self energy contribution of the level \a) 
due to interactions with level jcr'), corresponding to the 
diagram of Fig. [3] and given by 

S,,, (k, lUJn) ^~^Y.l 7^73^-' ^^n.) (8) 

X G[,?(q- k,ir2,„ - iw„). 

The large-A^ renormalized interaction vertex (the T- 
matrix) is determined by r~J, (q, ^$7™) = — — 

Caa'{^.,i^m), whcrc Caa' {^.^ i^n) IS the Correlator of 
the molecular (Cooper-pair) field operator B^rcr' (q) = 

/ I^ttT^ ^(t(p + q)V'cr'(~P)' ttiat to leading order in 1/A^ 
is given by 



C..,(q,zf]™) = ^ rfr e'"'"-(r, {b,,, (q, r)^^. (q, 0)}) 

= ^ E / G^"Hp + q, + *f^™)Gi']^(-p, -ZC.„) 

_ _ /■ d^V 1 - »F(^a,p+) - "-F(Cg',P-) 

I 



and p± = p±q/2, ^^,p = eCT,p-Ai<T, nF(a;) = l/(e'5^ -I- 1). 
An explicit expression for F^J, (q, uj) at zero temperature 
may be found in the Appendix of Ref. [39l . 

For our analysis, we will need the retarded fermionic 
self-energy S^^, (k, uj) at real frequencies and at fi- 
nite temperature. We denote with index R/A the re- 
tarded/advanced functions, i.e., functions analytical in 
the upper/lower half-planes of the complex frequency. In 
some cases, it can be obtained directly from Eo-o-' (k, ioJn) 
via a replacement iujn lo + i5. However, in general it 



is rather difficult to deal with discrete Matsubara sums. 
The approach we adopt here is to find the imaginary part 
of the retarded self-energy Im [S^^./ (k, w)] and obtain the 
real part via the Kramers-Kronig relation. 
Applying a Lehmann spectral representation 

f{iuj„) = - I dz ^— (10) 

to jl]), we find 



1 



1 



(11) 



After summing over the Matsubara frequencies Vtm and performing the standard analytical continuation iuj„ 
UJ -t- i5, we arrive to 



E«,,(k,c^) = -- 



d3q 



N J (27r)3 



dz 
2^ 



dz' nB{z) -\- np(z') 
27r z' — z + ll) + iS 



Im[F«,,(q,z)]lm G^/'"(q-k,/ 



(12) 



5 



where nsiz) = jin—i- The imaginary part of this function can be obtained by performing the z' integral to obtain 
(denoting z by Q) 



N J (27r)3 
Using the Kramers-Kronig relation 



^Im [r^,, (q, f])] Im [g^/") (q - k, - c^)] [ub + tif (f! - u;) 



Re [/(z)] = 



dz 



(13) 



(14) 



we obtain the real part of the self energy 



N J (27r)3 



2^ 



-Im[r^,,(q,17)]Re G^i'\ci-k,n 



+Re [r^,, (q, n)] Im [g;/°^ (q - k, 1] - c^)J nj.(f^ - uj) 
I 



(15) 



As mentioned in the Introduction, we are interested 
in the experimentally-relevant case where the number of 
atoms in each hyperfine state (species) is independently 
conserved. The Fermi wavevector /cpcr for each hyper- 
fine state is related to the corresponding density by the 
standard relation n^r = kp^/ (67r^). Since the atom densi- 
ties are fixed at ria-, the nontrivial self-energy corrections 
modify the dispersions and therefore lead to shifts in the 
chemical potentials fXa^ from their bare values of fia-o ■ 

More specifically, the chemical potentials are deter- 
mined by the condition G^^{k-pcr,^ = 0) = 0, yielding 
to lowest order in 1/N the relation 



5fia = Re [i;o.(fcFCT, 0)] , 



(16) 



where S^a are the chemical potential shifts for states a. 
At this point it is important to remark that although 
Eq. (fT6|) is strictly a self-consistent equation, where 
appears inside the self-energy, it would be overstepping 
the accuracy of the 1 /N expansion to use anything other 
than the bare values /io-o for the non-interacting Fermi 
gas at a particular polarization inside Sa-o-' • 



Another quantity of interest is the two-particle corre- 
lation function 



X,.'(q,*^^„) = - dr J d3re-'i-+-.Mr.{^t(r,^)^^,(r,r)V'I,(0,0)V,(0,0)}), 
where m„ = 27rn//3 is the bosonic Matsubara frequency. To order 1/N, we find for a =^ a' , 

P^J (27r)3 
1 1 .c-^ f d^'k' 

^^^'^zZ J {2-kY ^'^^^ ' iuj'()Ga' (k' -f q, ioj'i + iVn)Ta^a{k + k' + q, iujrn + iyj'i + ivn) 
which, upon using a Lehmann spectral representation and analytic continuation to real frequencies becomes 

A<,(k, L0)A„, (k + q, LO') ^ 



X«,,(q,z.)= / du I dw' 

5 



(27r)3 



v ~\- uj ~ uj' + iS 



1 

TV 



Wdzi 



(17) 



(18) 



1 



iwr^j g^^^C^^ ^^M. (k + q, z.)A.(k-, zM. (k- + q, z.) ; + ^' X 

nBF{z2,Zz,Z^) nBFizi,Z4,Z5) nB_F(zi,Z3,Z5) riBF{z2,Z4,Z5) 



Z2 — zi — V — i5 Zi — z^ — V — i5 



Z5 - Z2 - Z3 



Z5 - Zl - Z4 



Z5 — zi — Z3 — J/ — i(5 Z5 — Z2 — Z4 + J/ + i(5 

(19) 



6 



where nBF{x,y,z) 
np{x)np{ii), and 



[1 - npix) - npiu)] riBiz) 



A^(k,u;) = --lm[G^{k, 



(20) 



is the atomic spectral function. 

We will be particularly interested in the case where the 
state |3) is empty, i.e., fcpa = 0. In this case, it is useful to 
parametrize the two remaining Fermi wavevector fcp i and 
A:f2 in terms of the Fermi energy, ep, and the polarization, 
P. The Fermi energy is defined as ep = kp'^ /(2m), where 
n = rii + 712 = ^F'^/(37r^). The polarization, a measure 
of the population imbalance between states |1) and |2), 
is defined as 



P = 



"1 - "2 _ (fcpi)^ - {kF2f 



ni + n2 (fcpi)^ + (fcp2)^ 



(21) 



The evaluation of Eq. ([8]) is meaningful until the con- 
dition r~J, (q, 0) > is violated for some q, indicating 
an instability toward condensation of pairs with this mo- 
mentum. At unitarity the instability in ri2 occurs at po- 
larization P = 0.894(8), and q = 0.773. Although this is 
below the critical value Pc = 0.93 that follows from the 
mean field analysis, it should be remembered that the 
transition to the superfiuid state is first order so that the 
lower value corresponds to a spinodal instabilit y to some 
spatially inhomogeneous (FFLO) state [11, EE lH, El . 



III. FERMI LIQUID PROPERTIES: SPECTRAL 
FUNCTION 



One-particle quantities such as the momentum distri- 
butions Ua- (k) are conveniently expressed in terms of the 
fermion spectral functions, Ao-(k, cj), obtained from the 
retarded Green's function, Eq. [201 Note that in the case 
at hand, the medium and the interaction are isotropic 
and therefore the Green's function does not depend on 
the orientation of k. We will therefore drop the vector 
notation for k from now on. The momentum distribu- 
tion, n^{k) is then given by 



(k) — / duj np{uj)A„{k,ijj). 



(22) 



Near the Fermi surface, we can approximate the 
Green's function as a pole plus an incoherent spectrum 



G^{k,u;) = 



(23) 



W - i?a(fc) + «7cr(fc) 

where Ea-{k) is the renormalized spectrum of excitations, 
^a-{k) is the quasiparticle decay rate and Za- is the residue 
of the Green's function. This allows us to character- 
ize the corresponding strongly-interacting normal state 
in terms of the Fermi liquid parameters. The renormal- 
ized spectrum is connected to the real part of the self 
energy through the relation 

Ea{k) = e.fe ^fi + Re [l]^(fc, E^ik))] . (24) 



The quasiparticle decay rate is given by 70- (fc) 
— Im [S^ (fc, Ecr{k))] , whereas the residue is given by 



d_ 



Re[sf(fcp,,c^)] 



(25) 



The effective mass m* of the quasi particle at the Fermi 
surface is given by 



m 



k—k-pa- -t^— 



(26) 

We note that Eq. ([5]) coincides with the T-matrix ap- 
proximation, and therefore will reproduce the known re- 
sults in the dilute limit na^.^, ^ 1 obtained from the 
pseudopotential ^aaa'S{r — r'). 



1. Numerical Results 

For P = 0.9 at unitarity (i.e., 012 = 00), the results 
for the majority (|1)) and minority (|2)) momentum dis- 
tribution functions (Eq. are shown in Fig. [TJ The 
Migdal discontinuity at the Fermi wavevector is a mea- 
sure of the quasiparticle residue and correspondingly are 
given by Zi = 0.57 and Z2 — 0.30. An interesting fea- 
ture of this distribution is that ^2(0) ^ 0.45, showing 
strong depletion (from the maximum value of 1) of the 
single-particle states down to zero momentum. 

Quasiparticle properties, namely quasiparticle 
residues, effective masses and chemical potential shifts, 
are shown in Fig. [4] The properties are evaluated at large 
polarization P 0.9 — 1.0 and for densities appropriate 
to ^Li at the B — 838 G resonance. This corresponds 
to ai2 = 00, ai3 = — 3288ao and 023 = — 16080ao where 
flo is the Bohr radius {qq = 0.0529177nm). In order to 
make contact with the experiment of Schunck et al. , 
we set the atomic density such that fcpoia — —3.3 (which 
directly implies A;pa23 —16.1). However, we stress 
that, given that the hyperfine state |3) is empty, the 
single particles properties of the hyperfine states |1) and 
1 2) do not depend on the values of 023 and 013. In this 
sense their properties are universal and our numerical 
result are applicable to other fermionic resonant atoms, 
e.g., "^^K. On the other hand, the properties of the 
hyperfine state |3) are not universal in the sense that 
they depend on the interaction couplings. 

Using the formula Eq. (pS)) , we evaluate the quasiparti- 
cle residues Z^ at zero temperature. At full polarization 
(P = 1), the residue Zi = 1 is maximum since the state 
1 2) is empty. On the other hand, the residue Z2 = 0.47 
is much smaller than unity, indicating a large renormal- 
ization of atoms in state |2) due their interaction with 
atoms in state |1). 

The binding energy of the minority carrier is 112 = 
— 1.46ep, which can be compared to the Monte Carlo cal- 
culation of Ref. d^that found ^2 ~ — 0.93ep. 
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correlation by exciting atoms from occupied hyperfine 
states to another (usually empty) reference hyperfine 
state. As we will see below, the RF probe provides valu- 
able information about single particle excitations. 

In the RF experiments of Schunck et ai, one focuses on 
three different atomic hyperfine states of the ^Li atom. 
The two lowest states, |1) (majority) and |2) (minority), 
are populated and responsible for the superfluid corre- 
lations. The higher state, |3), is empty initially, and is 
used as a probe of atomic correlations in states |1) and 
2). An RF field, at sufficiently large frequency, is used 
to drive atoms from state |2) to state |3). 

In general the signal is expressed as a two-particle 
quantity. We calculate the rate for this process, de- 
scribed by a tunneling Hamiltonian that couples one of 
the species to a reference state |3) with a frequency lo 
detuned from the bare atomic transition. 

The RF field induces a transition between atomic levels 
primarily through the electronic spin. In the AC field 
of interest, the rotating wave approximation is used to 
described the tunneling Hamiltonian in terms of transfer 
matrix elements Vic.p between |2) to |3} state, i.e.. 



0.00 
-0.25 
-0.50 

£; -0.75 

=«■ -1.00 
-1.25 
-1.50 















___J2^_________-— 





0.90 0.92 0.94 0.96 0.98 1.00 



FIG. 4: Quasiparticle Fermi liquid properties: (A) residue Zcr 
(Eq. [25}, (B) effective mass m*(Eq. I26p and (C) the chemical 
potential shift Sjia (Eq. 116(1 of the hyperfine state |1) (dashed 
line), |2) (solid line) and |3) (dash-dotted line). 



(2^(2^ 



Vkpe-'^-'4{k)i;^{p) + h.c 



(27) 

For plane wave states, the tunneling matrix elements are 
Vkp = V{2Tr)^S{qL + k — p). Here ql ~ and ujl are the 
momentum and energy of the RF field. The conserva- 
tion of momentum in tunneling is to be contrasted with 
the analogous Hamiltonian used to model tunneling in 
superconducting-metal junction. 

The RF current is defined as /g = iVg = i[nT,N3] 
where N3 = J (fir ipl{r)ip^{r) is the number operator of 
the hyperfine state |3). Thus the current operator is 



hit) 



-iV 



(2^ 



^Uk)^^(k) - h.c. (28) 



It is important to note here that previous computa- 
tions 28] of the dispersion spectrum of the minority spin 
have all been done for the case of a single minority spin, 
in which case the problem is analogous to a Kondo/X- 
ray edge problem. Further, the dispersion spectrum is 
measured directly at k = 0. In contrast, we are consider- 
ing here the case of a finite density of minority spins, 
and the Fermi liquid properties near k = kF2, with 
\k — kp2\ ^ kp2- Even as kp2 0, this regime is distinct 
from the single minority spin case. With unitary inter- 
actions, it is not clear how the two regimes will connect. 
We leave this important issue open for future work. 

IV. RADIO FREQUENCY SPECTROSCOPY 

We turn next to the RF spectroscopy of Refs. H, [H. 
The RF spectroscopy is a technique used to probe atomic 



where as an accurate approximate we have neglected the 
extremely low RF photon momentum and set q^^ = 0. A 
standard linear response analysis gives 



(29) 



where v = uj^ + 012 — i^Js — ^J■2 + fJ-3 is the effective detun- 
ing. The retarded response function X^{0, v) is obtained 
from the analytic continuation of Eq. 1171 giving 



Im [X^(0, v)] = Im x'^y^^iQ, v 



Im 



X^\^^v) 



(30) 



where 
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Im 
Ini 



d^k f d^k' 



dbj [np{uj) — np{Lu + v)] Ao-(k, uj)Aa-' (k, u! + ly) 



(31) 



' — oo 



N J (27r)3 J (27r)3 

I — riBpii^ + iy,Lu' + iy,LU + Lu' + iy)Ga(k, LL!)Aa'(k, UJ + v)Ga{k' , Uj')Aa'{k! ., Uj' + v)I{a'a{k + k', + cj' + v) 

+ nBF(ijJ, ijj', + t^' + v)Ac{k, uj)Ga' (k, lj + i^)Aa{k' , u;')Ga' (k', + i')no-'o-(k + k' ,uj + u;' + 
+ 2 [npiuj) ~ npiuj + v)] Aa{k,uj)Acy>{k,uj + v)x 

[nsiio' +UJ + v)Gaik! , Lj')Go-'(k', t^' + v)W„i„{^ + k', a; + cj' + v) 
~ nF{oj')A,{k' , Lo')G,, (k', w' + v)T,,„{k + k!,LO + Lo' + v) 

- npiuj' + v)G„{k\ J)A„, (k', J + !/)r^'<,(k + k', + t^' + v)\ } (32) 

where n^rcr' (q, v) = -Tr^^Im [rg-g./ (q, v)\ . 

An explicit expression for r~J, (q, ^) at zero temperature may be found in the Appendix of Ref. [H. We note that 
in the hmit of a vanishing interaction between states |2) and |3), i.e., for 023 = 0, there are no vertex corrections to 
the A';^(0, v) correlator, and we find 



Im[X3^2(0,i^)] =lm[xi^'^(0,i.) 

d^k 



(27r)3 ]_ 



duj [nF{iL>) — npiijJ + v)] ^3(k, uj)A2ik, uo + v), 



(33) 



with all the interactions and corresponding 1/N corrections entering only through the atomic spectral functions 

^2,3 (k, u) . 

I 



2. Numerical Results 

We evaluate the RF spectrum for ^Li at the reso- 
nance from the above expressions, displaying the result in 
Fig. O As we mentioned in the previous section, interac- 
tions between particles introduce self-energy corrections 
to the fermions. In the calculation of the RF spectrum, 
it is important to consider all scattering lengths because 
the spectral weight of the states |2) and |3) are affected by 
all interactions. Following the formula Eq. ([3T|) . we can 
evaluate the contribution X^2\^i'^)- We find that the 
bosonic propagators ri3(q, w) and r23(q, w) have sharp 
poles at low momenta and positive energies. These pos- 
itive energy "Aleiner-Altshuler" poles are analogous to 
the bifermion mode discussed in the context of weak su- 
perconductors above the paramegnetic limit [44]. The 
sharp poles provide an additional spectral weight of the 
states |2), |3} but are not found to be significant contri- 
butions given the large energy gap associated with these 
modes. 



The interaction between states |2) and |3) introduces a 
vertex correction term, X^'^{0,i'), to the RF spectrum. 
Using the formula Eq. ([5^ and the condition that the 
reference state |3) is empty, we find that the vertex cor- 
rection is not significant in the limit P ^ 1. For exam- 
ple, for polarization P = 0.97 the vertex correction is less 
than 0.1 %, as shown in Fig. [S] 



In Figs. [H] and [21 we show the RF spectrum for the 
imbalanced polarization P = 0.93 and P = 0.97, re- 
spectively. Both RF signal peaks are positioned near 
A ^ 0.55 ep which is nearly equal to the chemical poten- 
tial difference between level |2) and |3) (see Fig. 2]). The 
results for the peak position compare favorably to the 
recent MIT experiment of Schunck et al. [l^ where they 
found A ~ 0.38 ep at finite temperature T/Tp ~ 0.08. 
The position of the peak for P = 0.93 is closer to the 
non-interacting peak than P — 0.97 case, and the in- 
creased pairing fluctuations at P = 0.93 further broaden 
the linewidth. 



V. CONCLUSION 

We have studied the normal state of a resonantly in- 
teracting, three component Fermi gas at large population 
imbalance. To this end, we utilized a well-controlled large 
N (atom flavor) expansion, with 1/A^ as the small expan- 
sion parameter and showed that it provides a systematic 
way to treat the strong interactions characteristic of the 
unitary scattering point of a Feshbach resonant system. 
Although the accuracy of the 1 /N expansion in the limit 
N ^ 1 remains unproven in general, it has provided rea- 
sonable estimates in other settings. Our main aim was 
to provide a qualitative interpretation of recent RF ex- 
periments [19], and to emphasize the importance of mea- 
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FIG. 5: RF spectrum for polarization P = 0.93: the intensity 
Isii^) (arbitrary units) vs. the detuning from the resonance 
frequency f measured in units of the Fermi energy ep. The 
peak frequency relative to non-interacting line (black line) is 
smaller than for the P — 0.97 case shown in Fig. [J] The 
linewidth broadens for lower polarizations due to the increase 
in pairing fluctuations. 
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FIG. 6: Vertex correction to the RF spectrum for P = 0.97. 
The axes are the same as in Fig. [5] The red line includes the 
vertex correction, while the black line does not. The contri- 
bution of the correction Xl^2'{0, u) is less than 0.1 %. 



suring quasiparticle properties to determine the nature 
of the ground state of imbalanced Fermi gases. Perform- 
ing the analysis to leading order in 1/N , we computed 
the atomic spectral functions and the momentum distri- 
bution functions with the characteristic Migdal discon- 
tinuity, thereby showing that the normal state is a con- 
ventional Fermi-liquid, albeit strongly renormalized. We 
then applied this formalism to analyze the RF excitation 
spectrum studied experimentally in Ref. [l^ Our con- 
clusion is that, indeed, the observed phenomenology can 
be well understood in terms of the strongly interacting, 
but conventional Fermi-liquid picture without resorting 
to exotic interpretations, such as, for example, pairing 
without condensation. 



Note added: Recent experiments by the MIT group 
[45I I have been interpreted in a physical picture consistent 
with that presented above, in contrast to their earlier 
interpretation of previous experiments . 
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